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CHAPTER  I 
INTRODUCTION 

The  object  of  this  study  is  to  analyze  the  high  frequency  radiated 
field  of  an  antenna  mounted  near  an  elliptic  disk.  The  elliptic  disk 
is  studied  because  of  its  generality  in  modelling  arbitrary  cross- 
section  flat  plates.  The  elliptic  disk  can  be  used  to  predict 
antenna  patterns  for  antennas  located  on  circular  or  elliptic  ground 
planes.  This  solution  can,  also,  be  used  to  model  an  endcap  of  a 
finite  elliptic  cylinder  used  for  modelling  structures  on  aircraft 
such  as  the  fuselage,  engines,  and  stores  or  structures  on  ships 
such  as  masks,  stores,  etc. 

The  basic  approach  taken  in  this  study  is  to  employ  the 
Geometrical  Theory  of  Diffraction  (GTD).  Since  GTD  is  a high  fre- 
quency method  the  lower  frequency  limit  of  this  solution  Is 
dictated  by  the  physical  dimensions  of  the  ellipse  and  the  source 
position.  The  radius  of  curvature  of  the  ellipse  should  be  greater 
than  a half-wavelength  and  the  source  should  be  at  least  a half- 
wavelength away  from  the  elliptic  edge. 

The  basic  problem  Is  analyzed  using  Keller's  Geometrical  Theory 
of  Diffraction  [1]  along  with  the  more  accurate  and  general  ray 
diffraction  coefficients  developed  at  The  Ohio  State  University  [2]. 

The  solution  for  the  elliptic  disk  case  is  modified  in  the  caustic 
and  pseudo-caustic  regions  by  employing  the  method  of  equivalent 
current  developed  by  Ryan  and  Peters  [3].  An  empirical  switching 
procedure  is  also  developed  in  order  to  determine  where  to  employ 
the  equivalent  current  solution  to  correct  these  caustics  to  yield 
a uniform  solution. 


CHAPTER  II 

THEORETICAL  BACKGROUND 

A.  INTRODUCTION 

The  radiated  field  of  an  antenna  positioned  near  an  elliptic 
disk  Is  analyzed  using  the  Geometrical  Theory  of  Diffraction  (GTD) . 
The  GTD  Is  a high  frequency  ray  optic  technique  which  allows  the 
radiated  field  to  be  analyzed  In  terms  of  individual  ray  components. 

In  this  study,  the  scattering  body  is  assumed  to  be  perfectly 
conducting  and  surrounded  entirely  by  free  space.  The  geometry 
associated  with  this  problem  is  Illustrated  in  Figure  1.  All  fields 
carry  an  exp(jut)  time  dependence  which  Is  assumed  and  suppressed. 
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B.  GEOMETRICAL  OPTICS  FIELDS 

The  incident  electric  field  E1  can  be  generated  by  an  arbitrary 
source.  In  this  study,  electric  or  magnetic  current  moments  pro- 
ducing a spherical  wave  are  assumed.  Note  that  an  arbitrary  source 
could  be  constructed  using  superposition  of  these  current  ele- 
ments [4].  The  far  field  pattern  of  the  source  is  written  as 


-jks 

E'(e,*)  = [e  F(e,*)  + « G(e,*)]  — 


(1) 


where  F(e,4>)  and  G( e ,<t> ) are  the  source  pattern  functions  and  s is 
the  distance  between  the  source  and  the  observation  point. 

The  reflected  electric  field  from  a perfectly  conducting  flat 
plate  is  determined  from  image  theory  and  is  expressed  as 


r(e,$)  = [e  Fr(e ,♦)  + $ Gr(e,*)] 


-jks' 


(2) 


where  s'  Is  the  distance  between  the  Image  source  and  the  ob- 
servation point.  The  geometry  associated  with  the  reflected  field 
is  shown  In  Figure  2. 

The  reflected  field  is  determined  using  the  source  field  in- 
cident on  the  point  of  reflection  Q and  the  boundary  conditions. 
The  source  field  incident  at  Q is  given  by 

Wl 


ti(ei,^1)  * [e  F(ei,^i)  + * G(e1  ,<(»i )] 


where  |I|  is  the  distance  between  the  source  and  the  point  of  re- 
flection and  the  angles  ei  and  ti  are  the  angles  of  incidence  from 
the  source  to  the  point  of  reflection.  The  boundary  conditions 
applied  at  Q are  given  by 


and 


r . i • r 1 

(3a) 

r • -t  • t 

(3b) 
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Figure  2.  Geometry  for  reflected  field 


where  z is  the  normal  to  the  surface  and  t is  a tangent  vector  to 
the  surface  defined  by 


Equations  (3a)  and  (3b)  can  be  rewritten  as 
(z  • e)Fr  + (z  • $)Gr  = (z  • 6j)F  + (z  • $j)G 


The  only  unknowns  in  the  above  equations  are  the  desired  functions 
Fr  and  Gr.  Solving  these  equations  yields 


pr  = (t  • j)hn  + (z  • »)At 


(z  • 0)(t  • 4)  - (z  • *)(t  • 0) 


(4a) 


Gr  = -C(z  • 9)Al  + (t  » »)An] 

(z  • 0)(t  • *)  - (z  • ♦)(€  • 0) 


where 

An  3 (z  • 0.)F  + (z  • 


A1  3 (t  • e^F  + (t  • ^.)G. 


C.  DIFFRACTION  BY  A WEDGE 

An  asymptotic  solution  for  the  diffraction  from  a conducting 
wedge  was  first  solved  by  Sommerfeld  [5].  Originally  plane  wave 
diffraction  coefficients  as  presented  by  Keller  [6]  were  used  as 
the  basis  for  the  GTD  solution;  however,  as  shown  in  Reference  [7], 
the  use  of  cylindrical  wave  diffraction  is  necessary  for  antenna 
applications.  Thus,  different  formulations  for  wedge  diffraction 
were  substituted  for  the  plane  wave  diffraction  coefficient  which 
is  the  basis  for  wedge  diffraction  theory.  Pauli  [8]  introduced 
the  Vg  function  as  a practical  formulation  for  a finite  angle  con- 
ducting wedge  Hutchins  and  Kouyoumjian  [9,10]  have  presented  a 
formula  ;jr  the  diffracted  field  which  significantly  improves  the 
accuracy  of  the  solution  over  that  obtained  from  Pauli's  form. 

This  improved  solution  provides  superior  results  in  the  transition 
regions  (near  the  incident  and  reflected  shadow  boundaries). 

The  three-dimensional  wedge  diffraction  geometry  is  depicted 
in  Figures  3 and  4.  The  source,  located  at  s'  (p ' ,<j>'  ,z' ),  generates 
a radiated  f field  given  by  F’(s).  The  source  can  be  an  arbitrary 
electric  or  magnetic  source  causing  cylindrical,  conical  or  spherical 
wave  incidence  on  the  wedge.  The  diffracted  vector  field  at  s(p,<j>,z) 
can  be  written  in  terms  of  a dyadic  diffraction  coefficient. 
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Figure  3.  The  basic  GTD  wedge  diffraction  problem 


Kouyoumjlan  and  Pathak  [11]  have  derived  a more  rigorous  basis  for 
the  GTD  formulation  and  have  shown  that  the  diffracted  field  may 
be  written  compactly  if  placed  in  terms  of  a ray  fixed  coordinate 
system.  This  ray  fixed  coordinate  system  is  centered  at  the  point 
of  diffraction  Q',  (or  points  of  diffraction  as  in  the  case  of 
plane  wave  incidence).  Q'  is  a unique  point  or  set  of  points  for  a 
given  source  and  observation  point.  The  incident  ray  diffracts  as 
a cone  of  rays  such  that  the  cone  half  angle  (e0)  equals  the  angle 
(6‘)  which  the  incident  ray  forms  with  the  edge  as  shown  in 
Figure  4. 

^Th§  orthoggnal  unit  vectors  associated  with  these  coordinates 
(s',B',$!,s,e_,$)  are  related  by 


where  I is  the  incident  direction  unit  vector  and  s is  the  unit 
vector  in  the  direction  of  diffraction. 
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Figure  4.  Geometry  for  three-dimensional  wedge 
diffraction  problem. 


The  diffracted  field  is  then  given  by 
^(s)  * E*(Q' ) • 0E(s,I)A(s)e'jks 


where 


°E  * -Ws  ' 


and 


"s 


-a 


2nj2nk  sin8„ 


+ cot 


F[kLa  ($-$')] 


'} 


<cot(I4n±  F[W-a+(W)] 

cot  ^^’))  F[kLa'(W)]JJ 


(5) 


where 


+ + 
a'(6)  * 2 cos2  j 


in  which  N is  the  Integer  which  comes  closest  to  satisfying  the 
following  equations 

2n*N  -8  * -w 


2nirN  -6  « -Hr 


with  a-  * V.  F(x)  is  defined  by 


8 


F(x)  = 2j|/n  ejxt  e 

>\*\ 


•jT2  dx 


and  1$  called  the  transition  function.  The  quantity  n is  defined  by 
the  wedge  angle  WA  where  WA  = (2-n)*.  For  the  elliptic  disk  case 
n = 2.  The  terms  are  associated  with  the  incident  field; 
whereas,  the  $+<(>'  terms  are  associated  with  the  reflected  field. 

The  quantities  A(s)  and  L are  defined  later. 

In  matrix  notation  the  diffracted  field  is  expressed  by  [9,10] 


Ed(s) 

-D  0 

E1 (Q  * ) 

II 

= 

s 

II 

Ed(s) 

_° 

e]_(Q') 

A(s)e 


-jks 


(6) 


The  Ds  coefficient  corresponds  to  the  it-field  component  parallel 
to  the  edge  with  the  boundary  condition  (acoustically  soft)  or 

(t|  ) - 0. 

'wedge 

The  Df,  coefficient  corresponds  to  the  it- field  component  perpendic- 
ular to  the  edge  with  the  boundary  (acoustically  hard)  or 


0. 


wedge J 


The  quantity  A(s)  in  Equation  (6)  is  a ray  divergent  factor 
given  in  general  by  [12] 


A(s) 


p 

sTp+sT 


i 


For  the  wedge  case  p*Pp  which  is  the  radius  of  curvature  of  the 
Incident  wavefront  In  the  plane  containing  the  incident  ray  and  the 
edge.  For  spherical  wave  incidence  A(s)  is  given  by 


A(s)  -J. 


s(s'  + s) 


L 


The  quantity  L Is  a distance  parameter  and  Is  given  In  general 
by  [12] 

2 

s(pc+s)p1p2  sin  e0 

L Pe(Pj+s)(p2+s) 

where  pj  and  pg  are  the  principal  radii  of  curvature  of  the  wavefront. 
For  a straight  edge,  L is  given  by  [12] 


s's  sinV 


s + s' 


(7) 


for  spherical  wave  Incidence. 

At  grazing  incidence  (♦'=0)  Ds  and  Dk  must  be  halved  since  the 
Incident  and  reflected  fields  merge  together  and  only  half  of  the 
total  field  on  the  surface  Is  the  incident  field  with  the  other 
half  being  the  reflected  field. 


CHAPTER  III 

SCATTERING  FROM  AN  ELLIPTIC  DISK  USING  THE 
GEOMETRICAL  THEORY  OF  DIFFRACTION 

A.  INTRODUCTION 

The  scattering  by  an  elliptic  disk  is  analyzed  in  this 
chapter  using  the  Geometrical  Theory  of  Diffraction  discussed  in 
Chapter  II.  The  geometry  of  the  disk  is  illustrated  in  Figure  1 
with  the  source  located  in  the  near  zone  of  the  disk,  and  the 
far  zone  scattered  fields  are  computed. 

B.  INCIDENT  FIELD 

The  incident  field  from  the  source  is  given  using  Equation  (1) 
by 

-jks 

E1  (e .♦)  = [0F(e ,$)  + ♦G((e .*)]  (8) 

^ A 

where  s ■ d + x$  • d. 


The  vector  d is  the  far  field  observation  direction  given  by 

d » cos  <t>  sin  6 x + sin  $ sin  0 y + cos  0 z 

and  djis  the  distance  from  the  origin  to  the  far  field  observer 
wltfli  xs  as  the  source  location  as  in  Figure  5.  The  phase  is 
referred  to  the  origin  using  the  usual  far  field  approximation 
yielding 

e-jks  -jkd  -jk  x • d 
— = — * *-3 — e 


It  is  possible  for  the  source  to  radiate  In  all  directions 
provided  that  the  ray  path  from  the  source  to  the  observation 
point  is  not  blocked  by  the  elliptic  disk.  If  the  source  ray  is 
blocked,  then  the  source  field  in  that  direction  is  set  equal  to 
zero  since  the  disk  is  assumed  to  be  perfectly  conducting. 
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Figure  5.  Intersection  of  incident  ray 
with  elliptic  disk. 


The  first  test  implemented  to  determine  if  the  observation 
direction  vector  (d)  and  the  vector  (xs)  from  the  center  of  the 
ellipse  to  the  source  are  on  opposite  sides  of  the  elliptic  disk 
as  ijj  Figure  5.  This  is  accomplished  by  comparing  the  sign  of 
z • xs  and  z • 3.  If  the  sign  of  these  two  dot  products  is  the 
same,  then  no  shadowing  can  occur  in  that  the  source  and  the  ob- 
server are  on  the  same  side  of  the  elliptic  disk.  If  the  sign 
Is  opposite,  further  testing  is  needed. 

For  the  next  test,  the  intersection  of  the  incident  ray  and 
the  plane  of  the  ellipse  is  found  by 


i 
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where  An  Is  the  orthogopal  distance  from  the  source  to, the  pljne  of 
the  ellipse  (An  » xs  • z)  and  Dn  Is  the  projection  of  d onto  z 
(Dn  » d • z). 

The  shadowing  test  Is  then  defined  such  that  If  1^1  > max(a,b) 
then  the  source  Is  not  shadowed.  If  |^|  < min(a,b)  shadowing  will 
occur.  If  mln(a.b)  < |x^ | < max(a,b)  a final  test  Is  implemented 
which  compares  |Ttl  with  the  exact  distance  from  the  center  of  the 
ellipse  to  the  edge  as  depicted  in  Figure  5.  The  angular  elliptic 
parameter  In  the  desired  direction  Is  found  using 


a \ * y 

b V * 


The  radial  distance  from  the  center  of  the  ellipse  to  the 
elliptic  edge  In  this  direction  is  found  from 

r~2  2 2 2 

Pe  s cos  vg  + b sin  vg  * 

Finally,  if  p?  < |x^ I no  shadowing  will  occur  and  if  pe  ^ Ix^l,  the 
Incident  ray  Ts  shadowed. 

C.  REFLECTED  FIELD 

Reflections  off  the  elliptic  disk  are  possible  only  if  the 
observation  point  and  the  source  are  located  on  the  same  side  of 
the  disk..  This  is  equivalent  to  saying  that  the  observation  di- 
rection (d)  and  the  vector  (xi)  from  the  center  of  the  ellipse  to 
the  image  location  are  on  opposite  sides  of  the  elliptic  disk  as 
in  Figure  2.  The  image  location  is  given  by 

aJ  a^  A 

X,  - xs  - 2 A„  2 


where  An  = xs  • z. 

The  procedure,  used  to  determine  if  a reflected  field  exists, 
parallels  the  section  on  source  shadowing  by  the  ellipse  except  for 
some  minor  changes.  The  reflected  field  exjsts  only  if  the  ray  from 
the  image  source  in  the  scatter  direction  (d)  intersects  the  elliptic 
disk.  If  this  intersection  does  not  occur  the  reflected  field  is 
set  equal  to  zero. 

If  it  is  determined  that  a reflected  field  exists,  then  the 
reflected  field  is  expressed  by 
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3 * 

*■•(•.♦>  ■ [®Fr(® ,♦)  ♦ «'(.,♦)]  e'Jk  Xi'd 

and  Fr(e .♦)  and  Gr(6,$)  can  be  found  using  Equations  (4a)  and  (4b). 
The  angles  of  incidence  8j  and  ^ from  the  source  to  the  point  of 
reflection  are  found  by  defining 


which  is  the  vector  frpm  the  source^ to  the  point  of  reflection 
located  at  aj  + An/Dn  d where  Dn  = d • z.  ei  and  are  then  found 
by  solving 


= tan 


D.  DIFFRACTED  FIELD 

The  diffracted  field  from  an  elliptic  disk  is  normally 
generated  from  discrete  points  along  the  elliptic  rim  which  satisfy 
the  law  of  diffraction.  This  law  states  that  at  the  diffraction 
point  (Q ‘ ) * the  angle  of  incidence  and  the  diffracted  cone  half- 
angle are  equal.  This  law  is  expressed  in  the  geometry  of  Figure 
6 as 


cos  8'q  = cos  eQ 


or 


I (v)  • e(v)  = d • e(v) 


(9) 


I (v)  is  the  unit  vector  from  the  source  to  the^ diffraction 
•).  If  Q*  is  located  at  the  point  (a  cos  v x ♦ b sin  v y) 
is  located  at  Ys  - xsx  + ysy  ♦ zsz  then 


where 
point  (Q 
and  the  source 
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LOCATION 


Figure  6.  Diffracted  field  geometry 


The  quantity  e(v)  Is  the  unit  vector  tangent  to  the  edge  at  Q 
is  given  by 


The  exact  solution  for  the  elliptic  parameter  v,  which  specifies 
the  location  of  the  points  of  diffraction,  is  complicated  for  an 
arbitrary  geometry.  The  solution,  in  general,  requires  solving  an 
eighth-order  polynomial  for  each  observation  direction  which  becomes 
very  time  consuming. 

Since  there  are,  in  general,  only  four  points  of  diffraction 
such  as  Illustrated  in  Figure  7,  a search  technique  has  been  used  to 


1 


Figure  7.  Curved  wedge  diffraction  points 
on  rim  of  elliptic  disk. 

find  these  points  of  diffraction  along  the  elliptic  rim.  l.<e  method 
is  based  upon  the  minimization  of  the  difference  between  the  two 
dot  products  given  in  Equation  (9).  An  error  term  is  given  by 

e(v)  « I(v)  • e(v)  - d • e(v). 


The  minimum  error  is  found  by  first  forming 

m e(Vj_-j)  + e(vj)  + e(vj+])  for  j = 0 to  360.  If 

|«j|  £ kj.il  and  <.  kj+il,  then  Vj  Is  the  closest  degree  to  a 

minimum  In  the  error.  In  order  for  this  minimum  to  be  a diffraction 
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point,  the  error  term  must  reduce  to  zero.  This  technique,  however, 
will  also  find  a minimum  in  the  error  which  is  not  a zero  crossing. 
In  order  to  remove  this  problem  and  to  determine  the  point  of  dif- 
fraction closer  than  one  degree  increments,  an  interpolation  scheme 
is  used  which  is  given  by 

vdiff  -V" 

where  Av  is  defined  as 

lv,  Itl  (e-d)  - (e-i) 

| e | + Ae  • I - |^j  (Ae*d)  - (e*I)  (e*d) 


and 


J g a - « 

Ae  * » -a  cos  v x - b sin  v z. 

If  |Av|  > 1 then  the  point  found  is  not  a diffraction  point  but  a 
minimum  In  the  error  function  and  is  ignored. 

The  location  of  these  diffraction  points  on  the  rim  of  the 
elliptic  disk  depends  upon  the  source  location,  the  observation 
direction,  and  the  physical  dimensions  of  the  ellipse.  If  the 
source  Is  located  on  the  z-axis  through  the  center  of  the  ellipse, 
as  in  Figure  8,  the  only  physical  dimensions  necessary  to  locate 


Figure  8.  Dipole  off  an  elliptic  disc. 
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the  diffraction  points  are  the  axial  ratio  and  the  height  ratio.  The 
axial  ratio  is  defined  as  the  ratio  a/b.  The  height  ratio  is 
defined  as  the  ratio  h/b.  The  quantities  a,  b,  and  h are  illustrated 
in  Figure  8. 

The  location  of  the  diffraction  points  on  an  ellipse  with  axial 
ratio  of  0.4  and  height  ratio  of  zero  for  a pattern  cut  in  the 
e ■ 90°  plane  (refer  to  Figure  8)  are  shown  in  Figure  9a.  In 
Figure  9b  four  observation  directions  were  chosen  from  Figure  9a 
to  demonstrate  where  the  points  of  diffraction  occur.  For  example 
at  4 = 25°  four  diffraction  points  occur  at  v * 10.5°,  135.6°, 

161.0,  and  233.9°  which  are  labeled  5,  6,  7 and  8,  respectively  on 
both  figures. 

The  location  of  the  points  of  diffraction  as  a function  of 
increasing  height  ratio,  i.e.,  as  the  antenna  is  moved  off  the 
disk  along  the  z-axis,  is  illustrated  In  Figure  10.  The  pattern  cut 
for  this  case  is  in  the  4=  90°  plane  for  an  ellipse  with  an  axial 
ratio  of  0.4. 


It  is  shown  in  Figure  11  how  the  location  of  these  diffraction 
points  are  effected  as  the  observation  direction  is  scanned  from  the 
4 = 90°  plane  to  the  ♦ = 180°  plane  in  10°  increments.  For  this  case 
the  axial  ratio  is  0.4  and  the  source  is  in  the  center  of  the  disk 
(height  ratio  is  zero).  Note  that  in  these  figures  if  0 and  ♦ are 
defined,  then  either  two  or  four  diffraction  points  exist. 

The  effect  of  changing  the  axial  ratio  is  shown  in  Figure  12 
for  an  ellipse  with  the  source  in  the  center  of  the  disk  and  the 
pattern  taken  in  the  ♦ = 90°  plane. 

After  the  discrete  points  of  diffraction  along  the  rim  have  been 
located,  the  diffracted  field  is  found  using  the  curved  wedge 
diffraction  coefficients  given  in  Chapter  II. 

The  diffracted  field  for  a single  diffraction  point  (Q‘)  on  the 
rim  of  the  elliptic  disk  illuminated  by  a current  moment  is  given 
using  Equation  (6)  by 


Ell 

s 

-D  0 
s 

En«n 

Fd 

LfiJ 

0 -D. 
L_  Q-j 

e[<q'2 

e-jkd 

3 


The  dyadic  diffraction  coefficients  for  this  case  are  found  using 
Equation  (5).  The  distance  parameter  L in  the  far  field  reduces 
to 
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Figure  9(a). 
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Location  of  diffraction  points  on  an  ellipse 
as  a function  of  observation  angle  (e=90°, 
axial  ratio  = 0.4,  height  ratio  » 0). 


0 (DEGREES) 


Figure  10.  Location  of  diffraction  points  as  a 
height  ratio  (4»=90° , axial  ratio  = 0 
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Figure  11(a).  Location  of  diffraction  points  as  a function  of 
observation  direction  (♦»90°-110°,  axial  ratio  « 
0.4,  height  ratio  ■ 0). 
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Figure  11(b).  Location  of  diffraction  points  as  a function  of 
observation  direction  (♦«110°-130°,  axial 
ratio  - 0.4,  height  ratio  ■ 0). 


Figure  11(c).  Location  of  diffraction  points  as  a function  of 
observation  direction  (♦■120°-150°,  axial  ratio 
■ 0.4,  height  ratio  ■ 0. 
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Figure  11(e).  Location  of  diffraction  points  as  a function  of 
observation  direction  (♦=170°-180°,  axial 
ratio  « 0.4,  height  ratio  * 0). 
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Figure  12(a).  Location  of  diffraction  points  as  a function  of 
axial  ratio  (♦“90°,  height  ratio  = 0). 


L * s'  sin  BQ 
where  s'  ■ |"f|. 


The  total  diffracted  field  is  obtained  by  summing  each  of  the 
diffracted  field  terms  provided  that  the  points  of  diffraction  are 
separated.  It  has  been  determined  that  if  the  diffraction  points  are 
separated  there  are  a maximum  of  four  such  points  which  will  satisfy 
the  laws  of  diffraction  such  as  shown  in  Figure  7.  This  four  point 
diffraction  solution  when  combined  with  the  incident  and  reflected 
fields  usually  predicts  patterns  which  are  in  very  good  agreement  with 
experimental  results.  However  for  certain  geometries  and  observation 
directions,  caustics  can  occur. 

The  axial  caustic  occurring  if  the  antenna  is  located  on  axis 
near  a circular  disk  (see  Figure  13)  has  been  treated  by  Ryan  and 


z 
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Figure  13.  Axial  caustics  of  a circular  disk. 

Peters  [13],  With  this  geometry  an  Infinite  number  of  diffraction 
points  exist  along  the  rim  of  the  disk  for  axial  observation, 
directions  (e  ■ 0 or  it).  The  infinite  number  of  diffraction  points 


implies  that  an  infinite  field  exists  in  the  neighborhood  of  e=  0 or 
n which  of  course,  is  not  the  true  physical  result.  However,  the 
equivalent  current  solution  given  in  the  next  section  can  be  used  to 
correct  the  diffracted  fields  in  these  caustic  regions  (i.e.,  the 
sectors  where  the  diffraction  solution  fails). 

A pseudo-caustic  problem  which  can  occur  with  the  elliptic  disk 
geometry  is  the  merging  of  diffraction  points  along  the  rim  as  the 
observation  direction  is  changed.  This  effect  is  illustrated  using 
the  configuration  shown  in  Figure  14a.  Solving  the  law  of  dif- 
fraction for  this  case  yields: 

A OBSERVATION 


I(v)  • e(v)  = d • e(v)  (10) 

where 

a a 

J _ a cos  v x + b sin  v y 

>1  a2cos2v  + b2  sin2v 

A A 

e = -a  sin  vx  + b cos  yy 
\|a2sin2v  + b2cos2v 

A A A 

d = sin  e x + cos  e z 

with  the  pattern  taken  in  the  y-z  plane  (♦=90°). 
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Equation  (10)  can  be  reduced  to  give 


b cos  v sin  e = 


2 2 

[b  -a  )sin  v cos  v 

? 2 2~ 
a cos  v + b sin  v 


If  cos  v * 0,  the  equality  holds;  therefore,  v = * w/2  are  two 
diffraction  points.  If  cos  v t 0 then  the  possibility  of  two  other 
diffraction  points  exists  and  Equation  (11)  becomes 


b sin  0 


[b2-a2)sln  v 

2 2 2 2~ 
a cos  v + b sin  v 


Squaring  both  sides  of  this  equation  and  solving  for  tan  v,  leads 
to 


tan  v 


+ ab  sin  o 


(b2-a2)2-b4sin2  e 


The  above  equation  has  two  real  solutions  for  v if 
(b2-a2)2  > b4  sin2  0 


which  implies 


| sin  0|  < |1  - (^)  | . 


If  the  inequality  is  true,  then  four  distinct  diffraction  points 
will  exist.  The  criterion  for  three  of  the  four  diffraction  points 
to  merge  into  one  occurs  if  a < /2  b for  this  pattern  cut. 

As  an  example  of  this  pseudo-caustic  effect,  the  case  where 
a a 2X  and  b * 5X  Is  chosen.  As  seen  Is  Figure  14b  caustics  occur 
at  the  angles  where  the  diffraction  points  merge  as  defined  by 
Equation  (12)  or 


t sin"1  {\  - {fyZy  t57.14°,  +122.86° . 
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Figure  14(b).  E-theta  radiation  pattern  in  the  y-z  plane 
for  the  configuration  in  Figure  14(a). 


The  merging  of  these  diffraction  points  as  well  as  their  position 
along  the  rim  is  seen  in  Figure  15.  If  the  pattern  is  taken  in 
the  x-z  plane,  which  is  similar  to  letting  a=5A  and  b=2x,  then 
a > /2  b and  no  caustics  occur  as  shown-in  Figure  14c. 


,Figure  14(c).  E-theta  radiation  pattern  for  the  x-z  plane 
for  the  configuration  in  Figure  14(a). 


The  plot  of  the  diffraction  points  for  a given  pattern  cut  is 
useful  in  showing  where  caustics  can  occur.  As  the  locus  of  dif- 
fraction points  approach  a vertical  line,  such  as  in  Figure  15,  a 
caustic  condition  is  being  approached. 
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Figure  15.  Diffraction  points  of  an  ellipse  as  a function  of 
observation  angle  (*=0°,  axial  ratio  = 0.4, 
height  ratio  = 0). 
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CHAPTER  IV 

CAUSTIC  CORRECTION  VIA  EQUIVALENT  CURRENT 


A.  INTRODUCTION 

The  equivalent  current  concept  was  implemented  by  Ryan  and 
Peters  [13]  to  correct  axial  caustics  occuring  in  the  back- 
scattering  of  axially  symmetric  cone  frustrums.  This  type  of 
caustic  occurs  because  GTD  predicts  infinite  fields  when  an 
infinite  number  of  diffracted  rays  occur.  The  term  equivalent  cur- 
rent is  used  because  the  currents  generated  are  a function  of  the 
scattering  direction  which  is  not  characteristic  of  true  line 
currents. 

The  equivalent  current  concept  is  used  to  correct  the  two 
types  of  problems  associated  with  the  GTD  four  point  diffraction 
solution  of  the  previous  section.  The  first  problem  is  the  axial 
type  caustic  similar  to  those  dealt  with  by  Ryan  and  Peters.  This 
case  occurs,  for  example,  with  a circular  disk  when  the  observation 
and  source  locations  are  both  axially  located.  The  second  problem 
is  a pseudo-caustic  effect  where  GTD  predicts  a discontinuity  in 
the  scattering  pattern  whenever  diffraction  points  merge  or 
disappear  from  view.  This  pseudo-caustic  effect  can  be  observed 
on  the  disk  when  three  of  the  four  diffraction  points  merge  into 
a single  diffraction  point  or  two  diffraction  points  merge  and 
disappear  as  the  observation  direction  is  changed  (see  Figure  15). 

B.  TECHNICAL  APPROACH 


The  equivalent  current  method  is  based  upon  the  two  dimensional 
diffraction  of  a plane  wave  incident  upon  a perfectly  conducting 
half  plane.  The  equivalent  current  is  found  by  comparing  the 
electric  and  magnetic  fields  generated  by  the  two  dimensional 
diffraction  process  at  each  point  along  the  edge  with  the  fields 
which  would  be  produced  by  equivalent  electric  and  magnetic 
current  elements.  The  diffracted  field  from  the  curved  edge  is 
found  by  numerically  integrating  these  equivalent  currents  in  the 
radiation  integral. 

Referring  to  Equation  (6),  the  two-dimensional  diffracted 
field  at  the  point  p is  given  by 


Ed(p) 


e!(Q')  d 


-jks 


VT 


(13a) 
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03b) 


Hd(p)  = H^(Q')  Dh 

V S 


In  terms  of  the  equivalent  electric  and  magnetic  l^ne  currents 
(Iz  and  M2)  along  the  edge,  these  fields  are 


‘d<P>  ■ -2o  *J  * e-j‘ 


H<I<P>  • ^ P’JkS  V 


Equating  these  two  expressions  implies  that 


E*(Q')  O(Q') 


(14a) 


«>)  °h(Q,)  |S7 


(14b) 


where 


^(Q’.6rt=90°) 


D _ ^h(Q'  »Bos90°)  Dfi 
h sin  60  ~ si% 


Ez(Q’)  and  Hz(Q')  are  the  electric  and  magnetic  components, 
respectively,  of  the  incident  field  parallel  to  the  edge. 


Using  the  above  equivalent  currents,  one  can  show  that  the 
diffracted  field  for  an  infinite  straight  edge  can  be  reproduced. 
This  is  done  using  the  following  line  integrals: 


1-d  _ 

4ir 


f -jks 

Jl(z')  sine  (z ' ) ^ — dz 


,d_  Jkzo  ell 

e z 


J¥  k<f>  S?  ^ 


3tt 

Ed  = ej  r' 
fce  e 


D|(z-)  dl,.  (15) 


The  above  integral  can  be  solved  using  a stationary  phase  argument 
since  a stationary  phase  point  occurs  in  the  integrand  of  Equation 
(15)  at  the  diffraction  point.  From  stationary  phase  [14], 


{' 1 f *’«"  f"Us> 


where 


f(z’)  = s + z'  cos$„ 


f'(z')  = -cos  e(z^)  + cosbq  = 0 


In  this  case. 


F(z* ) = Ej(z')  D|(z')  (z’>  1 ejkz’  cosBc 


wl  th 


s =/x2  + y2  + (z  - z*)2 


using  the  geometry  illustrated  in  Figure  16.  Therefore, 


. 3n 

d _ T 


E9  ’ * 


[E  Et(z,)p«  I eJk2'C0Se0 

\J7ir  zv  s'  s sTnT*  s 


. r~r~  ,-jks  -Jk2'coslso  -J  i 


I 


sin^S 


which  reduces  to 


z 5 D*  e 


-jks 


sin  eQ  s 


The  above  equation  checks  with  (13a)  in  that  e is  in  the  -z 
direction  and  secondly 


D_ . 


s1neQ  -% 


Returning  to  Equation  (15)  and  using  the  fact  that  in  the 
far  field: 


s * d - d • it 


it  follows  that 
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OBSERVATION  POINT 


Figure  16.  Geometry  for  three-dimensional  wedge 
diffraction  problem. 


9 


for  the  diffracted  field.  In  order  to  include  contributions  in 
the  integrand  for  directions  not  lying  on  the  cone  of  diffraction 
the  substitution  sin2g0  = sin  e0sin  e(z')  is  used  [15].  The 
exact  substitution  used  here  is  not  critical  to  the  answer  and 
this  substitution  was  chosen  to  enforce  reciprocity.  This  sub- 
stitution yields: 


For  the  elliptic  disk  case 


so  the  diffracted  field  becomes 


where  z'  is  the  tangent  to  the  elliptic  edge.  Using  duality,  the 
magnetic  field  becomes 


In  the  far  field 


n 


and 


? = zq  H x d 
which  implies 


_ rfi  a, 
zoH  'z 


*eikd‘*  . 

Also 

z/  - a *t  - -eJ>  ♦ eJ.J'  . 

Therefore 


Dj^  >A“dv  . 


The  D|  and  Dg  diffraction  coefficients  are  two-dimensional 
diffraction  coefficients  computed  by  assuming  a straight  edge 
which  Is  tangent  to  the  elliptic  edge  at  the  Integration  point. 
D|  and  Dg  are  computed  using  Equation  (5)  with  s1n$0  set  equal 
to  unity  which  implies  L ■ |T| . 
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The  above  equations  for  and  Ej  are  numerically  integrated 
using  the  following  technique: 


2n  f2l?  irhfvl 

f(x)dx  = M(x)eJ^x'dx  = 1 

o ■’o  x =0 

o 

in  steps 
Of  AX 

where  M(x)  and  $(*)  are  the  magnitude  and  phase  of  f(x),  re- 
spectively. Over  the  interval  Ax,  the  magnitude  is  assumed  constant 
and  the  phase  is  assumed  to  be  linear.  The  integral  now  becomes 


r2n  2ir 

| f (x)dx  % l f (x  ) Ax 

o x =0  0 

0 

in  steps 
of  Ax 


. AX 

f T 

AX 

T 


M(x)eJ*(x)dx 


where 


, dM*n. 


dx 


d>  ( xQ+Ax ) - ( xQ-Ax ) 
' 2Ax 


If  Ax  is  one  degree,  the  integral  can  be  reduced  to 


f 


f(x)dx  % 180 


180° 

I f(xj 


sin  A 


v°° 


o'  A 


where 


♦(x0+l°)  - 4(x0-l°) 


The  equivalent  current  solution  is  only  truly  valid  in  the 
vicinity  of  a caustic.  This  method  can,  however,  yield  good 
results  in  other  scatter  directions  due  to  the  stationary  phase 
nature  of  the  Integrand  at  points  of  diffraction.  One  region 
where  this  method  is  not  usually  valid  is  in  the  area  around  the 
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Incident  and  reflected  shadow  boundaries.  This  effect  can  be 
seen  by  examining  the  problem  Illustrated  In  Figure  17  and  also 
can  be  seen  as  a discontinuity  In  the  Integrand  of  the  equivalent 
current  solution  as  shown  In  Figure  18,  where  the  reflected  shadow 
boundary  Is  at  78.7°. 


OBSERVATION 


Figure  17(a).  Dipole  off  an  elliptic  disc. 


A computational  problem  associated  with  the  equivalent 
current  solution  Is  that  It  requires  Integration  around  the 
elliptic  edge  for  each  observation  direction.  This  requires  a 
large  amount  of  computer  time  when  compared  with  the  GTD  four  point 
diffraction  solution  which  requires  calculation  at  a maximum  of 
four  points  along  the  rim. 

Since  the  four  point  diffraction  solution  and  the  equivalent 
current  solution  fall  for  different  reasons,  the  discontinuities 
In  the  pattern  caused  by  the  diffracted  field  will.  In  general, 
appear  In  different  regions  of  the  pattern.  Thus,  the  two  methods 
may  be  combined,  in  most  cases,  to  yield  a smooth  pattern  which 
avoids  the  difficulties  associated  with  each  Individual  method. 

One  case  where  this  will  not  usually  yield  a smooth  pattern  Is 
when  the  caustic  and  a shadow  boundary  coincide. 
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Figure  17(b) 


Radiation  pattern  of  the  configuration  in 
Figure  17(b);  the  sharp  peaks  in  the  4-point 
solution  occur  at  the  pseudo-caustics  where 
the  equivalent  current  solution  must  be 
employed. 


INTEGRAND  MAGNITUDE  INTEGRAND  MAGNITUDE 


V OF  ELLIPSE 


Figure  18(a).  Dlscontlnutles  occurlng  In  the  Integrand 
magnitude  of  the  equivalent  current 
solution  near  the  shadow  boundary. 


The  four  point  diffraction  solution  is  used  whenever  possible 
since  it  is  the  more  numerically  efficient.  The  equivalent  current 
solution  is  used  to  correct  the  caustic  and  pseudo-caustic  problems 
associated  with  the  four  point  diffraction  solution. 

The  mechanism  chosen  to  switch  between  the  two  methods  is 
based  upon  the  stationary  phase  nature  of  the  integrand  in  the 
equivalent  current  solution  in  the  neighborhood  of  a diffraction 
point.  The  axial  caustic  problem,  associated  with  the  geometry 
illustrated  in  Figure  13,  is  apparent  from  the  phase  of  the 

integrand  as  seen  in  Figure  19.  Note  that  the  phase  is  constant 

at  the  caustics  (e=(f)  and  has  a broad  stationary  phase  peak  in 
the  vicinity  of  the  caustic.  The  phase  of  the  integrand  is  shown 
in  Figure  20  for  the  pseudo-caustic  case  of  Figure  14b.  Four 
distinct  phase  extremums  are  seen  when  e=45°  which  correspond 
to  the  four  diffraction  points  at  v=32°,  90°,  148°,  and  270°. 

As  the  observation  angle  changes  to  e*57°,  three  of  the  four  diffraction 
points  merge  into  one  at  v * 90°.  The  stationary  phase  peak  at  this 
point  becomes  very  broad  as  seen  in  Figure  20.  In  both  the  caustic 

and  pseudo-caustic  cases,  it  is  noted  that  the  width  of  these 

extremums  is  a good  indication  that  a caustic  is  present  and 
the  equivalent  current  solution  should  be  implemented. 

As  a result  of  the  above  arguments,  the  second  derivitive  of 
the  phase  at  the  diffraction  points  yields  a good  measure  of  the  width 
of  these  stationary  phase  peaks,  such  that  it  can  be  used  to  determine 
the  switching  procedure.  The  phase  of  the  integrand  can  be 
specified  in  a simplified  form  by  the  phase  of  the  diffracted 
field  at  the  diffraction  point  v^ . Thus 

"j  - -|f|  + d • R 

where 

A 

R = a cos  v^  x + b sin  v^  y 
and 

r.s-?s 


using  the  geometry  illustrated  in  Figure  21.  The  second 
derivitive  is  computed  numerically  by 


ni 


n1+1  ~ ni-l 
2av 
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Figure  21.  Geometry  for  switching  procedure. 


The  test  value  used  to  determine  when  the  equivalent  current 
solution  should  be  implemented,  was  found  by  empirical  methods,  i.e., 
deciding  where  the  switching  should  occur  for  several  test  patterns 
and  observing  the  value  of  ‘A*  . The  test  value  chosen  was  0.01 
(degrees )•! . Thus  if  |*ni|  is  less  than  0.01  (deg)-l  the  equivalent 
current  solution  is  used.  Otherwise,  the  four  point  diffraction 
solution  is  sufficient. 

An  application  of  this  switching  procedure  is  seen  using  the 
geometry  illustrated  in  Figure  17a.  The  pseudo-caustics  associated 
with  the  four  point  diffraction  solution  are  shown  in  Figure  17b. 

The  shadow  boundary  problems  associated  with  the  equivalent 
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current  solution  can  also  be  seen  in  Figure  17b.  The  combined 
solution  avoids  both  of  these  difficulties  by  employing  the 
equivalent  current  solution  only  in  the  pseudo-caustic  regions. 

The  combined  solution  is  shown  in  Figure  22. 

The  accuracy  of  the  combined  solution  can  be  demonstrated 
by  comparison  with  measured  results.  In  Figure  23a  the  axial 
caustics  associated  with  the  6TD  four  point  diffraction  solution 
is  observed  at  e=0°  and  180°.  The  combined  solution  which  cor- 
rects these  caustics  via  equivalent  current  is  seen  in  the  same 
figure  and  is  compared  with  measured  results  taken  on  a 3X  circular 
disk  [16]  as  shown  in  Figure  23b.  The  comparison  between  measured 
results  and  the  combined  solution  for  the  pseudo-caustic  case 
(see  Figures  14a  and  b)  can  be  seen  in  Figure  24. 


MEASURED  ON  CIRCULAR  Gnd  PLANE 
COMPOSITE  6TD  SOLUTION 


PATTERN  ANGLE  9 

Figure  23(b).  Comparison  of  combined  solution  with  measured 
results  on  a 3A  circular  disk. 


Figure  24(b).  E-theta  radiation  pattern  for  4>=45° 
corresponding  to  Figure  14(a). 


CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 

The  object  of  this  study  has  been  to  analyze  the  radiated 
field  of  an  antenna  mounted  near  an  elliptic  disk  at  high 
frequency.  The  basic  approach  used  in  this  study  to  analyze  the 
scattering  by  an  elliptic  disk  was  to  employ  the  Geometrical 
Theory  of  Diffraction  discussed  briefly  in  Chapter  II.  The 
incident,  reflected,  and  diffracted  fields  were  discussed  in 
Chapter  III.  The  diffracted  field  was  described  using  the  GTD 
four  point  diffraction  method. 

In  Chapter  IV  the  method  of  equivalent  currents  was  employed 
to  correct  for  caustics  and  pseudo-caustics  associated  with  the 
GTD  four  point  diffraction  solution.  An  empirical  switching 
procedure  was  introduced  in  order  to  apply  the  equivalent  current 
solution  whenever  caustics  occur  in  the  four  point  diffraction 
solution. 

The  combined  solution  for  the  elliptic  disk  was  compared  with 
measured  results  on  a circular  and  elliptic  disk  to  verify  its 
accuracy. 

The  equivalent  current  solution  has  been  shown  to  be  an 
effective  method  of  removing  caustics  associated  with  the  GTD 
four  point  diffraction  solution  on  the  elliptic  disk.  The 
combined  elliptic  disk  solution  using  the  GTD  four  point  dif- 
fraction solution  with  the  equivalent  current  solution  in  the 
caustic  regions,  can  be  used  effectively  to  predict  patterns  of 
sources  located  near  an  elliptic  disk. 

The  elliptic  disk  solution  is  useful  in  predicting  antenna 
patterns  for  antenna  located  on  circular  or  elliptic  gorund  planes. 

The  elliptic  disk  can  also  be  used  as  an  approximation  to  an 
arbitrary  cross -section  flat  plate  or  as  an  endcap  for  an  elliptic 
cylinder.  The  elliptic  cylinder  model  can  be  used  to  model 
structures  on  an  aircraft  such  as  the  fuselage,  engines,  and 
stores  or  structures  on  ships  such  as  masks,  stores,  etc. 

I M 
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